Resonant Pair Exchange and Percolation in the Disordered Hubbard Model 
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We show that the effect of disorder on a Mott-Hubbard insufator of interacting electrons produces 
a quantum phase transition from an antiferromagnetic to a paramagnetic phase governed by a per- 
colation of the singly occupied weakly disordered sites. Near the transition, we propose that a new 
type of defect, formed by the resonant exchange between a spin singlet and a doubly occupied site 
with an attractive disorder potential, plays an important role. These resonant pair exchange defects 
reduce the staggered magnetization but enhance the coupling of the two spins and produce charac- 
teristic signatures in the temperature dependent specific heat and the non-Curie spin susceptibility, 
at temperatures on the order the hopping t, higher than the typical exchange scale J. 



PACS numbers: 75.10.mv, 75.10.Nr, 71.27. +a 

Recent theoretical and experimental studies show that 
the interplay of random disorder with correlations in 
electronic systems can produce striking effects. Such 
studies have largely focused on the metal insulator 
transition^, 0, Q, where the problem still remains open 
from the theoretical perspective. Also of great interest is 
the question of magnetism in strongly correlated disor- 
dered systems. Besides the destruction of magnetic long 
range order with increasing disorder, we show in this Let- 
ter that disorder induces new types of magnetic coupling, 
absent in the ordered system. 

We address these questions by studying the effect 
of disorder on the magnetic properties of the Hubbard 
model. Besides its importance as a theoretical paradigm, 
this model is also relevant for the high T c superconduc- 
tors, where mobile carriers are produced by stoichiomet- 
ric doping that is also a source of quenched disorder. At 
half filling, in the absence of disorder, the model describes 
a Mott insulator with a spin at every site, and antifer- 
romagnetic (AF) coupling of order J ~ t 2 /U between 
spins at adjacent sites. We show that in the presence 
of disorder, it is also possible to have a novel type of 
magnetic coupling that we call resonant pair exchange 
(RPE), which operates between particular pairs of sites. 
This interaction is the outcome of the resonance to or- 
der t between two configurations of spins on adjacent 
sites: (a) a spin singlet formed by single spins on adja- 
cent sites and (b) a non-magnetic doublon formed by a 
doubly occupied and empty pair due to a binary potential 
close to ({7/2 and — {7/2) at two nearest neighbor sites. 
This sort of 'defect pair' occurs with a finite probability 
in a disordered system, and has several important char- 
acteristics: (i) Similar to two-level systems in glasses, 
these defects produce a characteristic maximum in the 
specific heat, (ii) The staggered spin susceptibility is 
suppressed because of the mixing of the singlet configura- 
tions with a non-magnetic configuration. However, inter- 



estingly, the non-Curie behavior persists to temperatures 
T ~ t, which is much higher than the kinetic exchange 
scale J ~ t 2 /U < t. (iii) The resonant tunneling pro- 
duces high kinetic energy on the bond connecting the two 
sites. This provides a source of noise in RPE defects that 
should be trackable in conductance noise experiments^]. 

These RPE defects are most active near the antiferro- 
magnetic (AF) to paramagnetic (PM) transition driven 
by increasing disorder. Strong disorder generates two 
types of sites: non-magnetic unoccupied or doubly occu- 
pied sites, and magnetic sites with a single spinp. We 
find that a percolation-based description then becomes 
possible, as electron hopping results in coupling between 
neighbouring magnetic sites; increasing the disorder leads 
to a decreasing number of magnetic sites, and eventually 
to a transition marking the loss of long range AF order. 
We supplement this with other calculational approaches, 
including an inhomogeneous Hartree Fock (HF) calcu- 
lation whose results are found to agree remarkably well 
with quantum Monte Carlo simulations[f|0]. 

Model: We start with the repulsive Hubbard model on 
a two-dimensional square lattice with site disorder: 



H = -t ( c L c j<? + c ]o c ™) + U^n l] n ii 



(1) 



where c ia .(ci a ) is the electron creation (annihilation) op- 
erator with spin a on site i and n ia = cT is the cor- 
responding density operator. The first term describes 
hopping of electrons between nearest neighbors with am- 
plitude t. The second term is the repulsive on-site in- 
teraction for doubly occupied sites with strength {7, fi is 
the chemical potential that determines the filling and Vi 
is a random potential chosen from a uniform distribution 
between — V and V. 
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The Large U limit: In the atomic limit t — 0, the 
competition between repulsive interactions and disorder 
produces site-dependent integral occupancies rii in the 
ground state. Sites with Vi > U/2 are unoccupied; and 
those with Vi < —U/2 are doubly occupied; neither of 
these two types of sites has a free spin. On the other 
hand, sites with \Vi\ < U/2 have rii = 1 with a free spin 
residing on each such site. Thus the fraction of randomly 
placed singly-occupied (magnetic) sites is x — U/2V, 
while the remaining sites are nonmagnetic. 

In the atomic limit, the spin degeneracy of the ground 
state is 2 xN , where N is the total number of sites. The 
effect of turning on a small hopping amplitude t is to lift 
this degeneracy and produce an AF coupling of magni- 
tude Jij between nearest-neighbor (nn) singly-occupied 
sites i and j, by the well-known mechanism of kinetic 
exchange. To second order in t, the coupling Jo = 4t 2 /U 
in the pure system is modified by the disorder difference 
SV — Vi — Vj on the two sites to 



Jo 



1-(8V) 2 /U 2 



(2) 



The effective leading-order Hamiltonian is then H = 
Jij Si ■ Sj where Jij is nonzero only if both sites 
are magnetic. 

Resonant Pair Exchange: It is evident from the above 
analysis that there will be some rare regions where the 
disorder at a pair of neighbouring magnetic sites Vi and 
Vj is such that the conditions \Vi — Vj — U\ < t <C U 
and \Vi/j \ ~ U/2 hold. In that case, the denominator in 
Eq.|2]becomes very large, and the perturbative expression 
is no longer valid. In fact, in this regime the electron 
hopping couples these pairs of sites to first order in t, and 
we show that it induces a new type of coupling, which 
we call resonant pair exchange (RPE). The RPE process 
differs qualitatively from normal kinetic exchange, and 
has important consequences for the thermodynamic and 
transport properties of the system. 

Consider states with 2 electrons on a pair of sites char- 
acterized by disorder parameters Vi, V 2 , when the hop- 
ping t — 0. Of the total of 6 states, there are 3 singlet 
(S = 0) states |l,l) s , |2,0), |0, 2). and 3 triplet (S = 1) 
states | |, T), | I, I), 1 1 , 1) t - Since the Hamiltonian con- 
serves total spin S, we examine each subspace separately. 
Of the three singlet states, two states |2,0) and |0, 2) in- 
volve unequal charges at each of the two sites, whereas 
one state |1, l)g = l/y/2(\ |,J.) — | Lf)) involves one 
electron on each site. We are interested in the case when 
one of the two unequal-charge states is nearly degenerate 
with |1, l)s- For specificity, let us take V\ = —V% = U/2. 
Then 5 of the six states (the 3 triplet states and 2 singlet 
states |1, l)s and |0, 2)) are degenerate with energy —U, 
while |2,0) has energy U. The primary effect of nonzero 
but small hopping is to mix the two degenerate singlets, 
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FIG. 1: Eigenstates for two sites, (a) Resonant pair ex- 
change: With Vi = — V2 = U/2; the difference in energy be- 
tween the RPE singlet ground state and triplet excited state 
is A ~ t. (b) Kinetic exchange: With Vi = Vi = 0; the energy 
difference between the magnetic singlet ground state and the 
triplet excited state is ~ J ~ t 2 /U. 



and produce the eigenstates 



m = 



1 

V2 



l,l)s±|0,2) 



(3) 



with eigenvalues —U — \[2t and — U + \ 2t respectively. 
The resulting pattern of energy levels is as shown in 
Fig. El 

In the case V\ — Vi =0, the states |2, 0) and |0, 2) have 
a high energy U , while the hopping lowers the energy of 
|1, 1)5 with respect to the triplet state by an energy of 
order t 2 /U (Fig. Qb). This is the mechanism of kinetic 
exchange, which produces AF correlations in the dimer 
ground state that persist up to T ~ t 2 /U. With reso- 
nant pair exchange, by contrast, the AF correlations in 
the ground state are reduced (as the non-magnetic state 
|0, 2) is mixed in); but these correlations persist to much 
higher T ~ t than in the case of kinetic exchange. As 
a consequence, both the uniform susceptibility x(l = 0) 
and the ordering susceptibility x(9 — 7r ) show deviations 
from Curie (free moment) behaviour for T < t, signal- 
ing the onset of AF correlations due to resonant pair 
exchange, as shown in Fig. Similar magnetic effects 
are expected whenever the conditions | U + Vi — Vj | < t 
and \Vi/j\ ~ U/2 are satisfied for neighboring pairs of 
sites (ij). The fraction of such pairs in the disordered 
Hubbard model is of the order of (t/2V) 2 . As the contri- 
bution to energy lowering from each pair is t, the overall 
contribution to the ground state energy is of order t 3 /V 2 , 
which is of higher order in t than from the majority of 
pairs, which are coupled by normal kinetic exchange. 
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FIG. 2: xM* an d x(0)t versus T/t for different U/t for two 
sites at half filling and with Vi = (7/2 and V 2 = -(7/2. The 
square symbol is for U = 8t and the circle symbol for U — 16t. 
When scaled by t, the curves for different T and U collapse 
to one curve. The inset depicts T\ versus \n(T/t), showing 
that x has Curie behavior at large T. 



Specific Heat and Spin Susceptibility: Since sites cou- 
pled by resonant pair exchange are relatively rare and 
have a very different level structure (see Fig.^l from the 
majority of pairs, they act as localized centers and give 
rise to a distinctive signature in the specific heat C v and 
Xavg, much as two level centers do in glasses. For an ar- 
bitrary pair, the energy splitting A between the singlet 
ground state and the triplet states is 



ht 2 



u -sv 



u -sv 



(4) 



where SV = \Vi — V2 1 - Here we have assumed the state 
with energy U+SV decouples from other two levels. If V\ 
and V2 are chosen from a uniform distribution between 
—V to V and we regard each pair as isolated from the 
others, the probability distribution for the splitting A is 
given by 



P(A) 



1 
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(5) 



The total specific heat from such pairs is C v — 
J^ mam dAP(A)c v (A) where A m i n and A rnax can be ob- 
tained by substituting and U for SV in Eq. 0] since 
the integration is only over singly occupied sites. Not- 
ing that the average energy for this two level system is 
E(A) = -Aexp(A/T)/[exp(A/T) + 3], the correspond- 
ing specific heat is c v (A) = dE/dT, shown in Fig. El 

We see that at T > t most of the contribution to the 
total specific heat is from RPE sites, due to their large 
splitting. In this regime C v varies as 1/T 2 . On lowering 




FIG. 3: The averaged susceptibility Xavg and specific heat 
versus the temperature. Both quantities show a peak at T ~ 
J where J ~ t 2 /U is the exchange coupling. The resonating 
pair exchange defects contribute at higher temperatures T ~ 
t. The distribution of the energy level splittings is shown in 
the inset and shows that it is dominated by small splittings 
~ J. 



T, C v has a peak at T = 2t 2 /XU with A ~ 2.85. Most 
of the peak weight comes from pairs with small splitting 
(A cx t 2 /U). In the low-T regime, C v decays exponen- 
tially as T — > 0; this form of the decay found within the 
pair approximation would change if the system supports 
extended (spin wave like) states. 

The above arguments hold also for the averaged sus- 
ceptibility. For the two level system mentioned above 
spin susceptibility is x + ~(A,T) = 2/T[exp(A/T) + 3]. 
The susceptibility averaged over pairs shows Curie be- 
havior (cx 1/T) at high T and a peak at T sa 2t 2 /U (see 
Fig. [SJ; below this temperature, triplet states make a 
very small contribution to the susceptibility. As T tends 
to zero xtvg ^ aus as ex P(~A/T). Therefore the high 
temperature (T > t) behavior of the specific heat and 
susceptibility is determined by the RPE sites while the 
low temperature behavior is governed by pairs of sites 
with small energy splittings. 

Percolation of Magnetic Sites: Recall that in the 
atomic limit, a fraction x — U/2V of sites is magnetic 
with rii = 1, while the remaining sites with m = or 2, 
are nonmagnetic. Since only nearest neighbour magnetic 
sites are coupled if t is small, the system maps to a dilute 
magnet @- 

A necessary condition for antifcrromagnetic long range 
order (AFLRO) is that there be an infinite connected 
cluster of singly occupied sites, i.e. x > x c where x c is the 
site percolation threshold for a given lattice geometry 0- 
Thus in the limit t/U <C 1, we can rigorously rule out 
LRO for x < x c . Further, in this limit, the model reduces 
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FIG. 4: Solid line is AF order parameter rrv vs 1 - i = 
1 — U /2V using HF approximation for a 2D square lattice 
of the model mentioned in Eq. Q the system size is 28 x 28, 
(7 = 4 and T = 0. The dashed and dotted lines are stag- 
gered magnetization vs dilution in the 2D quantum Heisen- 
berg model with spin- 1/2 using quantum Monte Carlo0, 
and spin wave theory 10] . The square symbols are experimen- 
tal data points using neutron scattering^. The percolation 
threshold is 1 — x c ~ 0.41. 

to a dilute Heisenberg AF with random couplings. The 
question is then whether there is AFLRO for x > x c . 
Since Jy in Eq. 2 is bounded below by Jo, and a dilute 
AF with uniform coupling between magnetic sites on a 
square lattice seems to show AFLRO |rj Ft Eof . we expect 
AFLRO in the disordered Hubbard model with t/U <C 1, 
provided that V < U/2x c . The occurrence of a small 
number of resonant pairs in the infinite cluster would 
result in a slight loss of AF order in the ground state. 
However, if the temperature is much larger than t 2 /U 
but still of the order of t, only the resonant pairs would 
make a non-Curie contribution to magnetic properties. 

Results of the Hartree-Fock approximation: The HF 
approximation, allowing for inhomogeneity in spin and 
densities on all the sites, sheds light on both magnetic 
and conducting properties of the disordered Hubbard 
model. Earlier work using this approximation focused 
on 3D ^lj and the effects of moving away from half fill- 
ing. Here we focus on the case of half filling in 2D. 
We find that adding site disorder leads to closing of 
the Mott gap in the density of states (DOS), while 
AFLRO persist upto the percolation threshold of mag- 
netic sites, 2V cr ~ £7/0.59 (Fig. BJ. At an intermediate 
regime of disorder and interaction the system has metal- 
lic behavior |l2j. Paramagnetic sites start appearing only 
when V > [7/2. For disorder strength V > V cr the sys- 
tem breaks into clusters of AF sites with no long range 
order, and displays a glassy behavior; the final state of 
self-consistent iterations depends on the initial inputs of 



the variational parameters of the trial HF Hamiltonian. 
With increasing disorder the size of AF clusters shrinks 
further, and at the limit of very large disorder the system 
is a paramagnetic Anderson insulator. 

Figure 0] shows a good consistency of the staggered 
magnetization obtained within the HF approximation 
with results of other theoretical and experimental stud- 
ies of the diluted Heisenberg model. The HF approxi- 
mation captures rather subtle effects as well. In partic- 
ular the next nn coupling, coming from a fourth order 
expansion in t can compete with the nn coupling to pro- 
duce occasional mismatches in the alignments of partic- 
ular spins/clusters — and these are reproduced by HF 
treatment. Further for a honeycomb lattice at half filling 
the AFLRO vanishes at a disorder strength, predicted by 
the percolation picture. 

Conclusion: There are two significant results in this 
paper: Firstly, we have identified a new type of disorder 
induced defect that involves pairs of binary potentials 
and produces a magnetic state with a reduced staggered 
magnetization but with an enhanced coupling. Secondly, 
we have shown that the disordered Hubbard model at 
large U can be mapped to a disordered diluted AF spin- 
1/2 quantum Heisenberg model. This allows us to con- 
nect the existence of AFLRO with a percolating infinite 
cluster of magnetic sites with \Vi\ < [7/2. We find a 
remarkable consistency between the percolation picture 
and HF results even for intermediate U. 
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